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FOREWORD 


This report describes the work performed by Scientific Systems, Inc. 
from 12/12/84 to 7/10/85 under National Aeronautics and Space 
Administration Award No. NAS1-17946. The work represents the first of a 
three-phase project aimed at development of theory, computational tech- 
niques, and a commercial grade control synthesis/analysis software package 
for large space structures; the approach is based on recent stable fac- 
torization and H«j techniques. The objective of Phase I effort was to 
demonstrate the feasibility of the concept. Dr. Earnest Armstrong was the 
project director for NASA. 

Dr. Hamid Razavi was the principal investigator. Professor M. 
Vidyasagar of the University of Waterloo was the consultant. Special 
thanks go to Ms. Alina Bernat for supervising the documentation. 



SUMMARY 


Large space structures are characterized by the features that (i) they 
are in general infinite-dimensional systems, and (ii) they have large num- 
bers of undamped or lightly damped poles. Any attempt to apply linear 
control theory to large space structures must therefore take into account 
these features. Phase I of this project consisted of an attempt to apply 
the recently developed Stable Factorization (SF) design philosophy to 
problems of large space structures, with particular attention to the 
aspects of robustness and fault tolerance. 

The final report on the Phase I effort consists of four sections, each 
devoted to one task. The first three sections report new theoretical 
results, while the last consists of a design example. Significant new 
results were obtained in all four tasks of the project. More specifically, 
(a) an innovative approach to order reduction was obtained, (b) stabilizing 
controller structures for plants with an infinite number of unstable poles 
were determined under some conditions, (c) conditions for simultaneous sta- 
bilizability of an infinite number of plants were explored, and (d) a 
fault-tolerant controller design that stabilizes a flexible structure model 
was obtained which is robust against one failure condition. 

The overall objective of Phase II effort will be to further develop the 
analytical foundation laid in Phase I, evolve computational algorithms, and 
test the stable factorization methodology on a realistic model of a large 
flexible structure. Thus by the end of Phase II effort the stage will be 
set for the development (in Phase III) of a state-of-the-art computer-aided 
control system synthesis/analysis package for large space structures, based 
on the stable factorization/H„ approach. 



0. INTRODUCTION 


In Phase I of this project, we proposed four tasks, namely (i) order 
reduction, (ii) coprime factorizations, (iii) robustness, and (iv) applica- 
tions. Of these, the first task was carried out to its natural completion, 
while the progress realized on the other three was more than sufficient to 
demonstrate their feasibility. 

We now describe in detail the four Phase I tasks and the results 
obtained. 

Throughout this report, the following symbols are used: S denotes the 

set of proper stable rational functions, i.e. functions that are proper and 
have all of their poles in the open left half-plane. A denotes the set of 
distributions f(.) of the form 

CO 

f(t) = I fi 6(t-t ± ) + f a (t), (0.1) 

i=0 

where 6 denotes the unit impulse distribution, 0<tg<ti< . . . . , f a ( ) is a 
Lebesgue measurable function, and in addition, 

CO CO 

Ilf 11 A = £ | fi j + j j f a (t) | dt (0.2) 

i=0 O' 

is finite. One can think of A as the set of all regular measures of 
bounded variation on the half-line [0,») that do not have any singular 
part; see Rudin (1962) or Hille and Phillips (1957). We refer to f a as the 
nonatomic part of f, and to the sum of delayed impulses as the atomic part 
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of f. The set A consists of all Laplace transforms £ of distributions in f, 
defined according to the familiar rule. 


f(s) = Z f ± e Sti + /“ f a (t) e st dt. (0.3) 

i=0 0 

If feA, we define IlfllA = Hf H a whete f is the inverse Laplace transform 
of f. Note that S is a subset of A, and that S, A, A are all rings. We 
use the symbols M(S), M(A) and M(A) to denote the set of matrices, of wha- 
tever order, whose elements all belong to S, A, A, respectively. 

On the set A it is sometimes convenient to use the so-called Hoo- norm , 
denoted by ll.lloo, and defined by 

Ilf lloo = sup j £(jm) | = sup | £(s) | . (0.4) 

cj Res>0 

Note that 

II f II oo < II £ HA for a ll fe A. (0.5) 

Finally, we introduce the so-called disk algebra , denoted by Ap, which 
consists of functions that are analytic on the open unit disk and con- 
tinuous on the closed unit disk. The norm on Ap is defined by 

Of II A h = max J f(z) | = max | f(eJ9)j . (0.6) 

I z I <1 0e[0 ,2ir] 
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1. TASK 1: ORDER REDUCTION 


Task 1 was as follows: find computational algorithms for the reduced- 

order modeling of large or infinite-dimensional systems, where the objec- 
tive is to minimize the graph metric distance between the true system and 
the reduced-order model. 

For ease of discussion, we reproduce below the definition and some pro- 
perties of the graph metric, as taken from Vidyasagar (1984), Vidyasagar 
(1985), and Vidyasagar and Kimura (1984). Given a matrix P, a pair (N,D) 
is said to be a right-coprime factorization (r.c.f.) of P if (i) N,DeM(A), 
(ii) P(s) = N(s)[D(s)]“l such that 

X(s) N(s) + Y(s) D(s) = I for all s. (i.i) 

If, in addition, it is true that 

N*(jco) N(joi) + D*(jw)D(ju>) = I for all u>, (1.2) 

where * denotes the conjugate transpose, then (N,D) is said to be 
normalized r.c.f. or n. r.c.f. of P. If P has an r.c.f., then it has an 
n. r.c.f.; see Sz. Nagy and Foias 1970. Note that (N,D) is an n. r.c.f. if 
and only if the matrix 



is inner . If P has an n. r.c.f., then it is unique to within right 
multiplication by an orthogonal matrix. 

Next, we define the graph metric. Let P 2 , P 2 be plants of the same 
dimensions, and let (N^jD^), (N 2 ,D 2 ) be n.r.c.f.'s of P]^ and P 2 , respec- 
tively. Define inner matrices B^,B 2 as in (3.1.3), and let 
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1IB 1 -B 2 Ulloo 


(1.4) 


6(Pi,P 2 ) 


inf 

UeM(A) , 


M U II 1 


Then the graph metric distance d(P^ ,P 2 ) between the two plants is defined 
by 


d(P 1 ,P 2 ) = max {fi(P 1 ,P 2 ), 6(P 2 ,P 1 )}. (1.5) 

Now we give some criteria for the robustness of feedback stability, 
given a plant P and a controller C of commensurate dimensions, define 


T(P,C) 


'(I+PC)~1 -P(I+CP)“1‘ 

C(I+PC)“1 -CP(I+CP)-1 



(I+PC)“1 [I -P]. 


( 1 . 6 ) 


Suppose a nominal plant-compensator pair (Pq.Cq) is stable, and suppose 
Pq,Co are perturbed to P,C respectively. Then the pair (P,C) is also 
stable provided 


d(P,P 0 ) HT(P 0 ,C 0 )lloo + d(C,C 0 ) IIT(C 0 ,P 0 ) lloo <1. (1.7) 

In particular, Cq stabilizes P provided 


d(P,P 0 ) II T( P 0 ,C 0 ) lloo <1. (1.8) 

This concludes the summary of known results, and also shows why the 
graph metric is a good measure of proximity in comparing unstable or 
nearly unstable systems. Note that there are no assumptions to the effect 
that P,Pq must have the same number of RHP poles, or that they should have 
no jco-axis poles, as is the case in the work of Doyle and Stein 1981. Thus 
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the above theory is ideally suited for studying the robustness of feedback 
stability in situations where the plant poles migrate across the stability 
boundary as the plant parameters vary. 

Now consider the problem of order reduction in the graph metric. 

Suppose P t is a "true" plant, for which we wish to construct an 
"approximate" model P a . Then P a is used to design a controller C, designed 
on the basis of the approximate model P a , to stabilize the true plant P t . 

We know this is the case provided 

d(P a ,P t ) < 1/ IT(P a , C) loo. (1.9) 

Next, suppose (Nj-.D^) is an n.r.c.f. for P^, and (N a ,D a ) is a (not 
necessarily normalized) r.c.f. of P a . Define 

( 1 . 10 ) 

V = II B t -M a II a, . (1.11) 

It is shown in Vidyasagar 1984 that if v<l, then 

d(P a ,P t ) < _2 v , (1.12) 

1-v 

In fact, a close examination of the proof reveals that (N a ,D a ) need only be 
a right fractional representation of P a ; it need not be coprime. This 
suggests the following procedure for approximating P t . 

(i) Find a normalized r.c.f. (Nt,Dt) for P t . 

(ii) Form the matrix Bf 

(iii) Approximate B t by M a . 



N t 



Na 

B t = 

_Pt_ 

9 

ii 

-Da. 
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To date, there are no known algorithms for optimal approximation in the 
norm II . Iloo. However, the pioneering work of Adamjan et al (1971, 1978) has 
led to a theory of optimal approximation in the so-called Hankel norm. 

More recently, Glover (1984) has derived bounds on the La,-norm of the 
approximation error B t -M a , where M a is found using the optimal Hankel norm 
technique. The method is described next: 

In the sequel, some assumptions are made concerning the plant Pt so 
that various technical assumptions are avoided. Task 2, discussed in the 
next section, is addressed to weakening some of these assumptions. Note 
that all of the assumptions below are automatically satisfied if P t is 
rational, so that the problem at hand is one of approximating a finite- 
dimensional plant by another of lower McMillan degree. 

(Al) The function Pt (s) has an unambiguous limit as | s | -*■» with Re 
s>0. Hence, we can speak about Pt(“>). 

(A2) Pt has only isolated singularities in the closed RHP, and each 
of these is a pole of finite McMillan degree. 

The above assumptions ensure (see Callier and Desoer 1978) that Pt has 
an n.r.c.f. B t . Moreover, since Pt(°°) is well-defined, so is Bt(°°)« Now 
introduce the bilinear transform 

z = s = i+£ . (1.13) 

s+1 1-z 

Then the function B t (s) gets transformed into an associated function 
B t [(l+z/(l-z) ] , which we again denote (by a slightly sloppy notation) as 
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B t (z). Note that when z=l, s=°°, so that B t (<») must be well-defined in 
order for this procedure to be valid. Now B t (z) belongs to the Hardy space 
Hoo over the unit disk (see Duren 1970). Hence it has the power series 

00 

B t (z) = Z B-jZ 1 . (1.14) 

i=0 

By the manner in which B t was formed, it is always a "tall” matrix. Let us 
square up B t as well as all of the Bj/s by adding columns of zeros, and 
form the block Hankel matrix 


'BO 

B 1 B 2 

e © • 



»1 

B 2 ... 




B 2 

• • • 


• 

(1.15) 

• 






n 

Now H can be viewed as a bounded operator on the space &2> consisting of 
all square summable, n-vector valued sequences. Moreover, H has finite 
rank if and only if B t is rational. Now consider the matrix H*H and exa- 
mine its spectrum. If Bt is rational (i.e. P t is lumped), then this 
spectrum consists solely of real eignevalues. In general, if Pt has an 
infinite-dimensional state space, the spectrum is real but need not contain 
only eigenvalues. Suppose we wish to approximate P t by a plant P a of 

McMillan degree k. If it happens that the spectrum of H*H consists of 
2 2 2 

eigenvalues aQ>a\> . . .> 0 -^, plus a subset of the real line contained in 
[0,a]J, this is easily done. Suppose first that k>n, and that 

2 2 2 2 2 2 
°0 >a l^ • • • ^ a k-n- >a k-n+l = a k-n+2 = ••• = °k. (1.16) 
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In this case, there exist n linearly independent eigenvectors x(^), x (^) t 
/ \ n 

x'- n -' , each belonging to £ 2 » °f the matrix H*H. That is, each vector 
x(l), ..., x( n ) has the form 


(i) (i) 

xC 1 ) = {xo , x 1 


U) 

x 0 

(i) 

X 1 


(1.17) 


Let y(i) denote Hx(i), and partition each y(i) as above. Finally, define 
the functions 


“ (i) 

Xi(z) = I x-jzi , 
i=0 


X(z) = [xjCz) .... x n (z)], (1.18) 

and note that X(z) is an nxn matrix-valued function. Define y±(z) and Y(z) 
accordingly. Then Adamjan et al. 1978 show that the function 


M a (z) = B t - o k X(z)[Y(z)]-l 

has McMillan degree k. Moreover, Glover 1984 has shown that 


(1.19) 


l!B t -M a !loo <21 Oi, (1.20) 

i=k+l 

if the spectrum of H*H consists solely of eigenvalues. If the multiplicity 
condition (1.16) does not hold, then one has to use so-called augmentation 
(see Adamjan et al. 1978). We will not discuss it here in the interest of 
brevity. 

So far we have accomplished the following: given a possibly distri- 

buted plant P t , we have formed the associated inner matrix B t , and found a 
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stable matrix M a of McMillan degree k that approximates B t . The last step 
is to go from the matrix M a to the approximating plant P a . This is done as 
follows: First, carry out the z to s transf ormation of (1.13), and par- 

tition M a (s) in the form 


M a (s) 


Then define 


N a (s) 

_ B a ( s )_ 


( 1 . 21 ) 


P a (s) = N a (s) [D a (s)]"l. (1.22) 

It can be shown that the McMillan degree of P a is at most k, and is generi- 
cally equal to k. An easy way to see this is as follows: Let 

(Aj^BjjijCjjijEia) be a state-space realization of the transfer matrix M a (s), 
and partition 



C N 


e n 

Cm — 


, E M = 



_ C D_ 

■ 

_ e d_ 


(1.23) 


Then, provided Ep is nonsingular, which it is generically, a state-space 
realization for P a is given by the quadruple 


A a = A^f-B^Ep 1 c D> B a = B M> 


(1.24) 


c a = c N -e N e D 1 c D> E a = E N E D 1 * 

In summary, we have presented an order reduction procedure that has the 
following novel and desirable features: 


(i) It allows approximation of distributed plants by lumped plants. 
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(ii) It allows one to specify a priori the McMillan degree of the 
approximate plant. 

(iii) It gives an a priori upper bound on the graph metric distance 
between the true and approximate plants. This bound, together 
with the stability condition (1.8), can be used to determine a 
priori whether a controller designed to stablize the approximate 
plant also stabilizes the true plant. 

One final feature of the approximation procedure is worth noting. Even 
if P t has infinitely many unstable poles, P a has only finitely many. 
Further, the locations of the unstable poles of P a need not correspond to 
those of P t . This feature defies conventional wisdom, which states that, 
when approximating an unstable plant, one must faithfully reproduce the 
"unstable part" of the plant. This is made possible by the use of the 
graph metric, which allows one to approximate the stable factors of a plant 
rather than the plant itself. 
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2. TASK 2: COPRIME FACTORIZATIONS 


Task 2 was as follows: Extend the existing theory of designing stabi- 

lizing controllers for linear distributed systems to the case where the 
plant has an infinite number of right half-plane poles, with particular 
emphasis on deriving conditions for the existence of a right-coprime fac- 
torization in such a case. 

The motivation behind this task is as follows: Vidyasagar 1978 showed 

that every rational matrix has both right and left coprime factorizations 
over A. More generally, Callier and Desoer 1978 introduced a class that 
they call fi , with the property that every matrix whose elements all belong 
to 6_ has both an r.c.f. and an l.c.f. However, functions in 8_ can have 
only finitely many singularities in the closed RHP, and each of these must 
be a pole of finite order. Nett et al 1983 showed that plants in M(fi_) 
could be stabilized by lumped compensators. The objective of this task is 
to see whether similar statements can be made about plants with infinitely 
many unstable poles (which, per force, cannot belong to &_) . 

To motivate the discussion, consider the plant 

p(s) = tanh s, (2.1) 

which has poles and zeros alternating along the jui-axis. Hence this is an 
infinite analogy of an LC admittance, i.e. a nondissipative system. It is 
easy to see that p(s) has a coprime factorization over A, since 
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p(s) = 


-2s 


( 2 . 2 ) 


1-e _ 

-2s 

1+e 

and since n(s)+d(s)=l. Hence it is certainly possible for plants with 
infinitely many unstable poles to have coprime factorizations. The objec- 
tive thus is to determine necessary and sufficient conditions for a plant p 
to have a coprime factorization. 

The question can be divided into two parts: (i) Given p(s), when does 
it have a fractional representation, i.e. when can it be expressed as a ratio 
n(s)/d(s), where n,dsA?; (ii) Given a p(s) with a fractional represen- 
tation, when does it have a coprime factorization? 

The first question is very difficult to answer in very general terms, 
since the number of possible functions p(s) is limitless. However, a few 
quick necessary conditions can be used to rule out functions that don * t 
have fractional representations. 

Theorem 2.1 A function p(s) can be expressed as a ratio n(s)/d(s) 
where n,deA only if it satisfies the following conditions: 

(i) In the open half-plane Re s>0, the only singularities of p are 
poles of finite order. 

(ii) On any vertical line o+ju), uieR, a>0, the zeros and poles of p are 
asymptotically almost periodic. 

Proof Any function feA has the form 
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00 


(2.3) 


f ( s) = 2 f± e st + f a (s), 

i=0 

where f a is the Laplace transform of an Lj function. Hence, by the 
Riemann-Lebesgue lemma, | f a (s) | ->-0 whenever | s j ->«> with Re s>0. Now the 
function 


f ap (ju>) = E f ± e~J a)t i (2.4) 

i=0 

is almost-periodic, as shown by Callier and Desoer 1972. In fact, if the 
delays t^ are periodic, i.e. if t^ = iT for some constant T>0 , then f ap is 
actually periodic. Hence, if p(s) = n(s)/d(s) and n,deA, then along any 
vertical line a+jto with a> 0, the functions n(a+jm) and d(o+jm) converge 
respectively to their almost periodic parts. This proves (ii). To prove 
(i), note that d is analytic in the open RHP, and as a consequence its 
zeros are isolated and of finite multiplicity. 

Using Theorem 2.1, one can quickly conclude that 

p (s) " (2,5) 

does not have a fractional representation, since its singularity at s=l is 
a branch point. 

The answer to the second question is, in principle, contained in the 
following result from Vidyasagar 1985, Section 8.1. 

Lemma 2.2 Suppose p(s) = n(s)/d(s), with n,deA. Then p has a coprime 
factorization if and only if the ideal in A generated by n and d is 
principal. 
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Though Lemma 2.2 contains a necessary as well as sufficient con- 
dition, it is difficult to apply it to a specific situation without further 
effort. The next result provides a useful sufficient condition. 

Theorem 2.3 Suppose p(s) has a fractional representation n(s)/d(s) 
where n,deA. Then p has a coprime factorization provided the number of 
RHP zeros of either n or d is finite. 


Proof We consider the case where n has only finitely many zeros in the 
closed RHP, the other case follows by symmetry. Let z^,...^^ denote the 
RHP zeros of n, and assume, by renumbering if necessary, that zi,.,.,z r are 
also zeros of d, while the rest are not. Let mj , . . . ,m r denote the minimum 
multiplicities of these as common zeros of n and d. Then, as shown by 
Callier and Desoer 1978, the function 


r 

f(s) = n 
i=l 

divides both n and d. Let 


|-S-Zi -t 


mi 


L s+i 

ni=n/m,d]=d/m. 


( 2 . 6 ) 

Then, by assumption, nj and d^ 


have no common zeros, and n^ is bounded away from zero at infinity; that 
is, if {s^jis any sequence in the closed RHP with j s^ | we have lim inf 
| n(s^) | >0. Hence 

inf | n^Cs) | + | di(s) j >0. (2.7) 

Re s>0 

This implies that nj and dj are coprime, as shown by Callier and Desoer 
1978. Since clearly p=n/d=nj7di, it follows that (n,d) is a coprime fac- 
torization of p. 
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We conclude this section by giving a very definitive theorem concerning 
coprime factorizations over the disk algebra Ap. 

Theorem 2.4 Suppose p(z) has a fractional representation n(z)/d(z), 
where n,deAp. Then p has a coprime factorization over Ap if and only if the 
following condition holds: Let Z n denote the set of zeros of n inside the 

closed unit disk, and define similarly. Then the cluster points of Z n 
and Z<j do not intersect. 

Before proceeding to the proof, let us explain what the theorem means. 
Recall that a point y is a cluster point of a set S if every neighborhood 
of y contains a point of S other than y. Also, note that since both n and 
d are analytic in the open unit disk, any zeros of these functions therein 
must be isolated. Hence the only cluster points of Z n and Z^ (if any) must 
be on the unit circle. Thus the theorem says that n and d can have common 
zeros in the open unit disk, and they can have common zeros on the unit 
circle provided that any such point is an isolated zero of at least one of 
the functions. But there cannot be a point on the unit circle which is a 
limit of both a sequence of zeros of n and a sequence of zeros of d. 

Proof "If" Suppose the hypothesis holds. Then there can only be fini- 
tely many common zeros of n and d in the closed unit disk, and each has 
only a finite multiplicity for at least one of the functions. Hence, as in 
the proof of Theorem 2.3, one can extract a polynomial as a common divi- 
sor of both functions. 

"Only if” This part of the proof closely follows Vidyasagar et al, 

1982. Suppose | zq j = 1, and that zq is the limit of a sequence {a-jJ of 


17 



zeros of n as well as of a sequence {b-jJ of zeros of d. We show that the 
ideal generated by n and d in Aq cannot be principal. Suppose by way of 
contradiction that there exist functions f,g,h,x,y in a D such that 


xn + yd = h, n = fh, d = gh. (2.8) 

Then gn=fd, so that f(ax)=0, g(b-^)=0. Also, since 

xf + yg = 1, (2.9) 

we have 

x(bi) f(bi) = 1 for all i (2.10) 

Now let i-*». Then lim bj=lim a^zg. Since f is continuous at zq, we have 

f ( zq) = lim f(a^) = 0. (2.11) 

i 

But (2.10) shows that 

x(zq) f (zq) = lim x(bi) f(b^) = 1. (2.12) 

This contradiction concludes the proof. 
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3. TASK 3: ROBUSTNESS 


Task 3 was as follows: Investigate the simultaneous stabilizability of 

an infinite family of systems. 

In its most abstract version, the problem is stated as follows in 
Vidyasagar 1985, Chapter 7. Let A be a first-countable topological space, 
let Aq be an element of A, and let P^, AeA be a family of plants. When 
does there exist a neighborhood B of Aq such that there exists a common 
stabilizing controller for all plants p^, AcB? The explanation behind this 
problem statement is as follows: A represents a set of physical parame- 

ters, and represents the corresponding system description. Aq is the 
nominal set of parameter values, and p^g is the corresponding nominal plant 
description. The objective is to determine whether a stabilizing 
controller for the nominal plant p^g also stabilizes all plants p^ when A is 
"sufficiently close” to Aq, i.e. whether there is a controller that achieves 
stability under nominal conditions and also maintains stability when the 
parameters are slightly perturbed. 

In one sense, the above problem formulation is incomplete. The objec- 
tive of controller design is not merely to stabilize the plant but also to 
achieve a desired closed-loop response. With this in mind, the following 
approach is used in Vidyasagar 1984, 1985. Define the closed-loop transfer 
matrix 
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H(P,C) 


'(I+PC) -1 

C(I+PC) _1 


-P(I+CP)~ r 

(I+CP) -1 


(3.1) 


Say that the family {p^} is robustly stabilizable at Ag if there exists a 
neighborhood B of Aq and a controller C such that 


(i) G stabilizes Px for all X in B, and 
(ii) H(Px ,C) is continuous in X at Xg. 


With this definition, a complete solution is given in Vidyasagar 1985, 
namely: robust stabilization is possible if and only if P^ is continuous 
in the graph topology at Xg. The graph topology is discussed at great 
length in this reference, and a characterization of continuity that is ade- 
quate for present purposes is this: P^ is continuous if and only if there 

exists an r.c.f. (Ng,Dg) of P\ and a family of r.c.f.’s (Nx>Dx) of Px such 
that N^Ng, Dx+Dg. 


In view of the aforementioned result, the only remaining interest in 
this problem concerns the case where the requirement that H(Px,C) be con- 
tinuous is removed. Suppose (Nx,Dx) is an r.c.f. of Px, that Nx,Dx are 
continuous, and approach Ng,Dg respectively as A+Ag. Let Pg=NgDg~l . Is 
there a controller that stabilizes Pg as well as Px for all X in some 
neighborhood of Xg? If Ng and Dg are right- coprime, the answer is clearly 
"yes”, by virtue of the aforecited results. Consider now the case where 
Ng »Dg are not coprime, so that there is pole-zero cancellation in the 
extended closed RHP. One can ask: Is it still possible to stabilize 

simultaneously Pg as well as Px? The next theorem shows that this is 
possible only if the cancellation is on the jw-axis or at infinity. 
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Theorem 3.1 Suppose N X ,D X are continuous with limits Nq,Dq, respec- 
tively as X-*-Xq. Suppose N X ,D X are right-coprime for X$Xq and let R be a 
greatest common right divisor of Nq,Dq. Let r = j R j , Pq = NqDq - !, P x = 
N X D X “1. If r has any zeros in the open RHP, then there does not exist a 
controller that stabilizes both Pq and P x , X+Xq. 

Proof We prove the contrapositive. Suppose C stabilizes both Pq and 
P x for X$Xq. Let Nq=NR, Dq=DR. Then N,D are right-coprime, and there 
exists a Jl.c.f. (Dq.Nq) of G such that 

D C D + N C N = I. (3.2) 

Now consider the return difference matrix 

U X = D C D X + N C N X . (3.3) 

Since C stabilizes P x , U x is unimodular for all X+Xo. Letting X>Xq in 
(3.3) gives 

R = lim U X . (3.4) 

X^-Xq 

In particular, r=j R | is the limit of | U X | . Now | U x | has all of its zeros 
in the open LHP, since it is a unit. Hence the zeros of r can at best be 
on the jw-axis or at infinity. 

As an application of Theorem 3.1, suppose 

P X ( s ) = l A - A , x 0 = 0, Po(s) = 1 (3.5) 

s-l+X 

Then no controller can stabilize both Pq and P x , X4=0. 
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Theorem 3.1 is the best possible result, because if { R j does not have 
open RHP zeros, then it may be possible to achieve simultaneous stabiliza- 
tion. This is illustrated by the following singular perturbations example. 

Let A = Ag=0, and consider the plants p\,Po described in state- 

space form by the equations 



y = x + z + 2u, (p x ) 


x = 3x-u, y = 2x+u, (pg) 


or in the transfer function form by 


p x (s) = 2 A s ^q -A) s - 5 f p o(s) = s-5 . 

Ax 2 +(l-A)s-3 s-3 

Now a coprime factorization of px is given by 


nx(s) = 2 As 2 + ( 1- A)s-5 > dx(s) = As 2 +(l-A)s-3 
( s+1 ) 2 (s+1) 2 


As A->0 , we get 


n A’ >n 0 = — — > d A' >d 0 = — - 

(s+1) 2 (s+1) 2 


Clearly ng and dg are not coprime, since they have the nominal common 
divisor 


r(s) = . J , . 
s+1 

However, the controller 

c(s) = - 8 

s+12 


stabilizes pg as well as p\ for sufficiently small positive values of A. 
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In this case the only zero of the common factor r is at infinity. 
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4. TASK 4: FLEXIBLE STRUCTURE EXAMPLE 


In order to illustrate the usefulness of the Stable Factorization 
approach in a specific control example, a four-disk model of a flexible struc- 
ture will be presented here. 




Aetuaeor 

Figure 1 

The disks are connected by torsion springs. Let denote the spring 
constant for the spring between i-th and (i+l)-th disk; let denote the 

inertia for i-th disk. The equations of motion for the system are: 
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where 0 £ denotes the angle of deflection of the i-th disk from normal posi- 
tion. Each dot indicates differentiation with respect to time. The right 
hand side indicates the scaled torques applied (in this case to the 3rd 
disk) . 

The motivation of this model is that by applying actuator inputs 
(torques) u-^ and disk i, one wishes to control the outputs 9 j , for i equal 
to or different from j. The transfer function from each input u^ to each 
output 0j can be easily calculated from the dynamic system (S). For the 
purpose of illustration we assume that the spring constants Ki are equal to 
unity. The nominal plant is further assumed to have Ji = 1, i=l, 2, 3, 4 
(e.g. equal inertia). Note that due to the fact that there is no damping 
in the system, the transfer function from any input to any output will have 
all of its poles and zeros located on the jar-axis. 

For example, the transfer function from u 3 to 64 will be of the form 
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PCs) = 


b(s) 
a (s) 


S 2 + g 2 
4 

S 2 ir (S 2 + a 2 i) • 
i=l 

where 0^ and aq are real positive scalars, corresponding to a zero at 
Z 1 = j3l and poles p^ = joii, i = 1, 2, 3, 4. The rigid body modes 
correspond to the poles s 2 = 0. 

In the following we present a fixed controller that simultaneously sta- 
bilizes a nominal plant Pj and a contingent plant P 2 . This is an example 
of the simultaneous stabilization techniques introduced in Vidyasagar and 
Viswanadham (1982) as discussed earlier. The particular pole zero loca- 
tions correspond to the experimental set-up described in Cannon and 
Rosenthal (1984), also reported in Razavi, Mehra, and Vidyasagar (1985). 

The computation of the actual controller using state-space versions of the 
stable factorization approach is carried out based on methods of Minto 
1985. 

The system under study consists of four steel disks connected by three 
steel rods as shown in Figure 1 . A control input is applied to the second 
disk from the top, and the two outputs are the positions of the first and 
fourth disks, respectively. Two different situations are considered. In 
the first, all four disks are identical. In the second, the first three 
disks are identical and the last disk has a moment of inertia that is one 


- 26 



quarter of those of the rest. In the nominal case, the system is denoted 
by Pi(s) and equals 

(4.1) 


(4.2) 

When the last disk is smaller than the rest, the system is 

b2i(s) 

p 2 ( s ) = 827,200 , (4.3) 

a 2 (s) |_b 2 2 (s)_ 

where 

a 2 (s) = S 2 (s±j29 . 34) (s±j21.48) (s+jll.938), 
b 2i(s) = bn(s), (4.4) 

b 22 (s) = (s±j25 .95) . 

The objective is to find a single controller C that stabilizes both pj 
and p 2 . In the present instance both systems are entirely undamped, i.e. 
all poles and zeros are purely imaginary. If we stabilize the two systems 
by moving the poles into the left half-plane, it is possible to come up 
with designs that yield closed-loop systems that are nominally stable, but 
have poles with very little damping. To rule out this possibility, we 


bn(s) 

Pl (s) = 827,200 

a l(s) _b 12 (s) 


ai ( s ) = s 2 (s±j 23.81) (s±j 18.22) (s±j9.863), 
bn(s) = (s±j 20.85) (s±j7 .96) 


bnCsl = ( s±i 1 2.8861 
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define the stability region to be the region { s : Re s<-l}. Thus, in order 

for the closed-loop system to be stable in this more restrictive sense, its 
poles must have real parts less than -1. 

Using the methods of Minto 1985, we have designed a seventeenth order 
controller that achieves the above objectives. The closed-loop poles of 
the two compensated plants are shown in Table 1. The controller itself is 
shown in modal form in Appendix 1. Note that modes 7,8,11, and 15 can 
essentially be ignored, resulting in a thirteenth order controller. 



(Pi) 


(P2> 


- 1 . 00110+00 
-9. 99910-01 
-1. 001 iD+0© 
-9. 99910-01 
-1. 00790+00 
-1. 00790+00 

-a. ooooo+oi 

-9, 99920-01 
-9. 99920-01 
-9. *4700-01 
-l. 73030+00 
-l. 73030+00 
-2. 20430900 

-a. *20*0+00 

—4. 04740+0© 
-4. 9*490+00 
-9.993*0+00 
-7. 00990+00 
-7. 99710+00 
-1. 03070+01 
-l. 01090+01 
-1. 01090+01 
-9. 734*0+00 
- 1 . 00020+01 
-9. 73480+00 


+2. 22940+0 li 
+2. 0*900+01 i 
-2. 22940+011 
-2. 08900+0 li 
+ 1. 2*370+011 
-1. 2*370+011 
♦2. 17440— I6i 
+7. 94390+001 
-7. 94390+001 
+?. 18810-101 
+«. 44070-01 1 
-8. 44070-011 
-3. 34230-091 
+2. 30920-091 
+2. 21140-101 
-7. *9990-101 
+9. 39430-101 
-l. 90770-101 
+3. 28390-111 
-4. laoto-ioi 
+2. 93230-011 
-2. 93230-01 1 
+ 2 . 10020-011 
+1. 47190-101 
- 2 . 10020-011 


-4. 03950-01 
-1 00000+01 
-1. 76160+00 
-l. 76160+00 
-1. 50150+00 
-2. 49750+00 
-3. 498*0+00 
-4. 50340+00 
-9. 49740+00 
>4. 90070+00 
-7. 49990+00 
-8. 90000+00 
-2 94940+01 
-2. 94940+01 
-2. 90170+0© 
-2. 90170+00 
- 1 . 00110+00 
-9. 99910-01 
- 1 . 00110+00 
-9. 99910-01 
-1. 00790+00 
-1. 00790+00 
-9. 9992D-01 
-9. 99920-01 
-4. 03990-01 


♦■3 9562D+00 1 
+ 1. 29970-121 
-3. 59410-01 1 
♦a. 39410-011 
+6. 35600- lOi 
-l. 3710D-091 
-3. 28450-091 
+4. 42460-091 
-2. 23430-091 
+4. 82830-101 
-3. 92340-11 i 
♦4. 49470-131 

-I. 19190+011 
+1. 19190+011 
-2. 79490+011 
+2. 79490+01 1 
-2. 22940+01 i 
-2. 08900+011 
+2. 22940+011 
+2. 08900+011 
-l. 28370+011 
+ 1. 28370+011 
-7 96390+001 
+7. 94390+001 
-9. 9 9420+00 i 


Table 1. Closed loop poles of nominal Cpi ) and 
contingent (P 2 ) plants 
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APPENDIX 1 


Controller Equations for Simultaneous Stabilization 

(Section 4) 


Controller dynamic equations: 

x c (t) = (AC)x c (t)+(BC)u c (t) 
dt 

y c ( t) = (CC)x c H) + (EC)u c (t) 


Controller input u c : 
output y c : 
state x c : 


2x1 vector 
lxl scalar 
17x1 vector 


Controller matrices: 


AC: 17x17 matrix (diagonal) 

BC: 17x2 matrix 

CC: 1x17 vector 

EC: 1x2 vector 


AC = DC - I 17 

1^7 : identity matrix of dimension 17 


A-l 



Diagonal elements of DC = 


-3. 16240+01 *1. 94990*01 i 
“3. 162*0*01 -1. 94990*011 
-1. 99870*00 *2. 77920*011 
-1. 99870*00 -2* 77920*011 
1.44090*01 *1. 18290*01 1 
1.44090*01 — l. 18290*011 
9.20470-01 *7.84940*001 
9.20470-01 -7.84940*001 
-1.84410*01 *4.44880*001 
-1.84410*01 -4.44880*001 
-7. 94470*00 *3. 34700*001 
-7. 94470*00 -3. 34700*001 
-9. 33100*00 *1. 97940*001 
-9. 33100*00 -1. 97940*001 
-1.00000*01 *3.47340-141 
- 8 . 99210*00 *4. 12080-011 
-8.99210*00 -4. 12080-011 


BC = 


-4. 78940-02 -t. 02200-021 
-1.09290*01 *2.17100*011 

1.70110*01 *4.62490*001 
-1.23110-03 *7. 19970-041 
-7. 97820-09 *7. 19990-041 
-1.18490*01 -3.26380*011 

2. 12770*01 *6. 21400*011 

2. 19890-09 -2. 28100-091 
- 8 . 42120-02 *2. 40380-021 
-1. 04490-01 -6. 99440-021 

4. 78220-02 *7. 23440-021 
-9. 81920-09 *3. 70440-071 

3. 97270*00 *9. 99470*001 
-7.48940-01 -4.83710-011 

9.41730*01 *1.64820*021 

3.41190*00 -1.27990*001 
-3. 99840-02 -3. 79310-021 


A-2 



cc = 


CaiuMAt 1 te 3 

-3. 97*30+01 -1. 38070+011 

Colgans 4 to 6 

73370+03 -4. 73330+031 

CoIimnio 7 to 9 

1- 30300+01 **9. 27330+001 

CoIimnio 10 to 13 
-7. 31980+00 +3.44890+011 

CoIimnio 13 to 19 
4.44070+01 -3. 88440+01 1 

CoIimnio 14 to 17 
-t. 13770-01 +3.00410-011 


-1. 73100+03 +3. 32940+031 
2. 48060+00 +1. 34390+011 

-2.00740+01 +1.93910+011 
-1.07140+04 -1.34130+031 
-4. 29870+01 +4. 99480+011 
2.34940+00 -7.99000-011 


-9. 32290+02 +8. 78980+031 
-1. 40430+01 +9. 83300+001 

3. 29370+01 -4. 23770+011 
3. 17330+00 +3. 10910+001 
-1.83490-03 -3. 43090-011 


EC = 


4. 47380-01 0. OOOCO+OO 
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